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MULTIVARIATE RBR TOOL AGING ADJUSTER 

RELATED PATENT APPLICATION 

U.S. Patent Application Serial # 10/313,501 (TSMC-01-1509), to the same inventors, 
filed on Dec. 6, 2002. 

BACKGROUND OF THE INVENTION 

(1) FIELD OF THE INVENTION 

The present invention relates to the fabrication of integrated circuit devices, and more 
particularly, to a method of monitoring tool performance using a multivariate advanced process 
control system in the fabrication of integrated circuit devices. 

(2) DESCRIPTION OF THE PRIOR ART 

Run by Run (RBR) feedback control schemes have been used widely to control and 
monitor the semiconductor manufacturing processes. Statistical process control is combined with 
feedback control and uses data from past runs to adjust the recipe for the next run. This method 
offers the potential for reducing process variability caused by equipment aging, chemical depletion, 
or fluctuation in ambient conditions with a minimal cost. Most literature in the RBR control, 
however, is only applicable to the manufacturing process with a single product by a single tool. 
For a foundry industry, due to different customer requirements, one of the most popular production 
styles is multiple products by multiple tools. For example, if a process runs eight products by five 
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tools, then there are forty combinations which will lead to a very complex situation using 
conventional RBR control schemes. If separate feedback schemes are set up for each combination, 
besides affecting the efficiency of the RBR control schemes, this method also ignores the useful 
information between the tools and products. For example, under a specific tool, an individual 
control scheme does feedback only on a specific product and it loses the opportunity to adjust for 
other products produced by the same tool. In addition, since customer orders may not be 
continuous and scheduled tool set-up may change fi*om time to time, the feedback scheme for any 
production lot cannot be successfully applied to the next production lot. Therefore, it is desired to 
provide a multiple run by run control process. 

A number of patents address process control issues. For example, U.S. Patent 5,963,902 
to Wang xises an expectation maximization (EM) model in a pattern recognition process. U.S. 
Patent 6,238,936 to Yu teaches mapping critical dimension (CD) of predetermined features at 
various process stages. A user can look at the data to see if adjustments should be made to an 
etching process. U.S. Patent 5,757,673 to Osheiski et al teaches automated analysis of data from 
self-tests of steppers uploaded as an ASCII file. U.S. Patent 6,245,581 to Bonser et al teaches a 
method of measuring critical dimension, analyzing the measurement, and performing a secondary 
etch process based on the analysis. U.S. Patent 5,862,054 to Li shows a process parameter 
monitoring system. Statistical Process Control (SPC) analysis is performed on data firom each 
machine. Results are displayed in graphical form. U.S. Patent 6,263,255 to Tan et al discloses a 
software framework for measurement, processing, and strategy control of a factory. U.S. Patent 
5,243,377 to Umatate et al shows a lithography control system. 
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SUMMARY OF THE INVENTION 



Accordingly, it is a primary object of the invention to provide a multiple run by run 
process control method. 

Another object of the invention is to provide a multiple run by run process control method 
that works across tools and products. 

Yet another object is to provide a multiple run by run process to detect tool aging. 

Yet another object is to provide a multiple run by run process to detect photolithography 
tool aging. 

A further object is to provide a multiple run by run process to detect photolithography tool 
aging based on exposure dose. 

A still further object is to provide a multiple run by run process to update a recipe based 
on past data. 

Another object of the invention is to provide a multiple run by run process to update a 
recipe for a particular tool and product based on tool decay and variation shift of product. 
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In accordance with the objects of the invention, a multiple run by run process is achieved. 
A plurality of tools and a plurality of products to be run on the tools are provided. A desired recipe 
is calculated for each product on each tool based on a previous recipe used for each product on each 
tool and deviation of a parameter from a target value. Thereafter, the plurality of products is run 
on the plurality of tools. The desired recipe is updated after each run of each tool. 

BRIEF DESCRIPTION OF THE DRAWINGS 

In the accompanying drawings forming a material part of this description, there is shown: 

Figs. 2 and 3 are cross-sectional representations of a photolithographic process. 

Fig. 3 is a diagram of the two phases of the advanced process control system of the 
present invention. 

Fig. 4 illustrates a flow chart of the feedback control portion of the statistical process 
control portion of Phase 2 of the advanced process control system of the present invention. 

Fig. 5 illustrates a flow chart shovsdng the gradual and rapid mode corrections of the 
advanced process control system of the present invention. 



4 



TS0M510 



DESCRIPTION OF THE PREFERRED EMBODIMENTS 

The process of the present invention is a muhiple run by run advanced process control 
system. This system is a model-based process control system which combines the advantages of 
statistical process control (SPC) with feed backward control to monitor processes and prescribe 
control action as needed. The APC system of the present invention uses statistical methodology to 
find the process optimization model and apply the process recipe to maintain the best process 
performance during a process shift or drift of controllable scope. The APC system diagnoses 
process variation and judges whether or not corrective action is to be taken. The process recipe is 
updated when necessary to meet the process target. Statistical methods are applied to calculate 
process variation to achieve a stable and accurate feedback algorithm. The process of the invention 
improves the process capability index by reducing the rework rate, increasing throughput, reducing 
consumption of materials (such as photoresist in a photolithography process), and decreasing 
engineer labor. 

The multiple run by run (MRBR) APC method of the present invention can be used for 
any process in semiconductor manufacturing. The photolithography process will be used to 
demonstrate the MRBR scheme of the present invention. It will be understood by those skilled in 
the art that the MRBR scheme of the present invention can be used in any process in semiconductor 
manufacturing and that the invention should not be limited to the photolithography process and the 
energy parameter described as an example herein. 
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Figs. 1 and 2 illustrate the photolithography process flow. Fig. 1 shows a wafer 10 
covered by a photoresist layer 20. A mask 30 is placed over the wafer. The photoresist is exposed 
to light 40 through the mask. The photoresist is developed to leave the patterned photoresist 
layer shown in Fig. 2. 

The quality characteristic of the photolithography process is critical dimension (CD), 
shown in Fig. 2. According to physical knowledge, CD is a function of the input recipe and the 
exposure dose or energy (E). The increment or decrement of the energy will lead to a 
proportional change in CD. Based on this relationship, the MRBR controller is to adjust the energy 
setting such that CD will meet a specific target. 

Fig. 3 is a diagram of the two phases of the MRBR system of the present invention. In a 
first phase 31, data is collected and entered into a database. This phase takes place off-line. First 
we build a prediction model for the photolithography process as follows: 

^run - (Eb/R + EReticIe) X Tooloflfset + ^correction (1) 

where Erun = process energy 
Eb/r= basic record energy 
EReticIe = mask error energy 
Tooloffset = tool energy ratio 
Ccorrection = k shift energy correction 
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The parameters above are defined for each process, reticle, and tool. Other information to 
be collected includes CD target, process time, and so on. Now, two datasets are to be constructed: 
Desired Energy and Adjusted Energy. These datasets are two-dimensional arrays D = [dy] and A = 
[ay], respectively, where i = l,...,n represents reticle ID and j = l,...,m represents tool ID. The input 
data contain several variables including mask level, tool ID, product ID, reticle ID, target, CD 
output, and energy. The adjusted energy is to compensate for the bias between output CD and 
target CD. By linear approximation, the following formula can be used to obtain the adjusted 
energy of each production run: 

Adjusted Energy = slope x (CD - Target) (2) 

The constant slope is assigned by the user, depending upon mask level, the magnitude of the final 
energy, and so on. Hence, the desired energy of each production run is given by: 

Desired Energy = Final Energy + Adjusted Energy (3) 
The value dy is equal to the original final energy plus the adjusted energy; that is: 

dy = ey H- ay, for i-l,...,n; j=l,...,m (4) 

Now the two tables are to be generated. For example, this procedure can be as follows: 
1. Input the values of the slope. 
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2. Data sorting by tool ID, mask level, target, and reticle ID. 

3. Use equations (1) and (2) to obtain the Adjusted Energy and Desired Energy of each production 

run. 

4. Desired Energy Table: 

a. Average the Desired Energy in the same tool ID and reticle ID 

b. Decompose the average Desired Energy dataset by tool ID into several sub-datasets. 

c. Merge the decomposed sub-datasets by reticle ID and set the value of missing cells to zero. 

5. Adjusted Energy Table: 

a. Average the Adjusted Energy in the same tool ID and reticle ID 

b. Decompose the average Adjusted Energy dataset by tool ID into several sub-datasets. 

c. Merge the decomposed sub-datasets by reticle ID and set the value of missing cells to zero. 

Now, we need to compensate for missing data in the two tables. We can use an 
expectation maximization (EM) algorithm, a recursive algorithm used to forecast or estimate 
the Desired Energy of zero-containing cells in the dataset. 

We may use the following procedure, for example: 

1 . Create a dummy array U = [uij] having the same size as the Desired Energy array. Let uy = 0 if 
dy is not missing (>0) and uy = 1 if dy is missing (=0). 

2. Calculate the missing numbers in each row and each column: 

m 

Ui. = Zuy, i=l,2,...,n (5) 
j=l 
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Uj = EuijJ=l,2,...m (6) 



3. Calculate the row means of the array D: 



X dij/(n-Uj) 



4. Replace every missing cell of the array D by the corresponding row means. 

5. Set the stop condition of the recursive algorithm; for example. Stop =105 and iteration = 30. 



6. If the stop condition is false, do: 

a. Calculate colunm means and row means: 



n 

d.j[k] = Zdij[k]/n,j=l,...,m 
i=l 



m 

di.[k] = Z djj[k]/m,i=l,...n 
j=l 



where k is the iteration number 



b. from i=l to n 
fromj = 1 to m 

Ifuij = l,then 

dij[k+l] = (m X (n X dj[k] - dij[k]) + n x (m x di.[k] - dij[k]) - (n x m x d„[k] - 
dij[k]))/((n-l)(m-l)) 
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c. Let M =X (dij[k+l] - dij[k])2, for all uy = 1. 

d. If M > Stop or k < iteration, go to a. 



For example, the desired energy table is given in Table 1 A. Two cells are empty. 



Table lA 

Product Tool 1 Tool 2 Tool 3 Tool 4 

1 31.38 28.24 33.17 34.96 

2 0 28.81 33.46 33.65 

3 29.11 0 31.51 32.72 



step 3: Calculate row means: 

(31.38+28.24+33.17+34.96)74 = 31.94 
(28.81+33.46+33.65)73 = 31.97 
(29. 11+31 .5 1+32.72)/3 = 31.11 



step 4: Replace missing cells with row means: 

Table IB 

Product Tool 1 Tool 2 Tool 3 Tool 4 

1 31.38 28.24 33.17 34.96 

2 31.97 28.81 33.46 33.65 

3 29.11 31.11 31.51 32.72 

step 6: Iterate. 
k=l: 

Row means: 
(3 1 .38+28.24+33 . 1 7+34.96)74 = 3 1 .94 
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(31.97+28.81+33.46+33.65)/4 = 31.97 

(29.1 1+31.1 1+31.51+32.72)/4 = 31.11 
Column means: 

(31.38+31.97+29.1 1)/3 = 30.82 

(28.24+28.81+31.1 1)/3 = 29.39 

(33.17+33.46+31.51)/3 = 32.71 

(34.96+33. 65+32.72)/3 = 33.78 
d21[2] = (4 X (3 X 30.82 - 31.97) + 3 x (4 x 31.97 - 31.97) - (12 x 31.67 - 31.97))/6 = 30.27 
d32[2] = (4x (3 X 29.39 -31.11) + 3 X (4x31.11 -31.11)- (12 x 31.67 - 31.1 1))/6 == 26.82 



Table IC 



Product 



Tool 1 
31.38 
30.27 
29.11 



Tool 2 
28.24 
28.81 
26.82 



Tool 3 
33.17 
33.46 
31.51 



Tool 4 
34.96 
33.65 
32.72 



2 
3 



k=2; Repeat 



Table ID 



Product 



Tool 1 
31.38 
30.98 
29.11 



Tool 2 
28.24 
28.81 
26.70 



Tool 3 
33.17 
33.46 
31.51 



Tool 4 
34.96 
33.65 
32.72 



2 
3 
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k=3; Repeat 



Table IE 



Product 



1 
2 
3 



Tool 1 
31.38 
31.00 
29.11 



Tool 2 
28.24 
28.81 
26.70 



Tool 3 
33.17 
33.46 
31.51 



Tool 4 
34.96 
33.65 
32.72 



Now, we need to obtain the optimal setting for each tool ratio using the reticle energy and 
the correct energy table. The correct energy table is generated by finding the difference between 
the desired energy and the actual energy for each tool and product combination. 

The ideal model for energy is 



This is equation (1) without the correction factor. The optimal energy setting could be modularized 
by tool ratio, base energy and reticle energy. In practice, there is a difference between the energy 
of the model and the optimal energy; this is the correction factor in equation (1). Ignoring the 

correction factor for now, we take the logarithm of both sides of equation (7) and get: 

log(Erun) =l0g(EB/R + EReticle) + log(Tooloffset) (8) 

To estimate the tool ratios or offsets, the following procedure may be used: 
1 . Take the logarithm, Idy = log(dij) 



Erun - (Eb/R+ EReticle) X Tooloffset 



(7) 
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2. Estimate the tool ratio: 

a. Suppose that the subscript index of the golden Tool ID is index j = jO. (The golden Tool is 
that tool which has the most stable quality in the process.) 

Let Fold be the original tool ratio of the golden tool. 

n 

temp-ratio = X (Idijo - leijo) x (l-Uijo)/(n-Ujo) 
1=1 

Then, Fnew =r old x exp(temp-ratio) 



b. Calculate tool ratios FjJ = l,...,m 

n 

Calculate the column means: Id.j = X Idi j In 

i=l 



Calculate tool ratios except j = jO: 
Fj = exp(ldj - Id.jo) x F new 
3. Calculate reticle energy: 

a. Input base energy Eb/r. 

m 

b. Calculate row means: Idj = ( X Idy) / m, i=l,...,n 

j=l 

m 

c. Let tti = Idi. - (Z Id.j )/m + Id.jo - log(r „ew) - log(EB/R) 

j=l 



Then, the Energy of the Reticle is 
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Ei = (exp(ai) - 1) X Eb/r, where Eb/r is the corresponding Base Energy of the Reticle 

ID. 

d. Now, create an array for recording the correct energy table, using the correction factor. 
C = [cij], i=l,...,n;j=l,...,m 

Cy = dij - (Eb/R + Ereticle) X Tooloffset 

This completes the first phase. 

The second phase of the APC system of the present invention is performed on-line. This 
phase includes adjusting the recipe by an Exponential Weighted Moving Average (EWMA) 
controller 34 and detecting the changes on tool effect and product effect using integrated statistical 
process control 33. The feed backward system of the present invention is shown in the flowchart of 
Fig. 4. 

First as shown in step 41, energy data is collected. This is the actual exposure energy 
used in the process step. Now, in step 43, calculate the Desired Energy where 

Desired Energy (New Final Energy) = Final Energy + w * slope * (Target - CD) 

That is, Eij(t) = Eij(t-1) + w x slope x (Target - CD) 
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The final energy is equal to the desired energy in Phase 1. The initial value was found in 
Phase 1. 

Now, in step 44, tool ratios are calculated. For each control ID, a tool aging ratio is 
calculated using the equation: 

Qt(i) = Qt.i(i)+log(Desired Energy)-log(Final Energy) 

where Q is the tool aging ratio and where Qo(i) = 0. 

In step 46, the tool ratio is checked against user-defined aging criteria. If the tool ratio 
does not meet the criteria, the tool ratio is re-calculated in step 48. If tool aging is detected, the tool 
ratio or tool offset is recalculated. 

In step 52, if no tool aging is detected, the tool ratio is not recalculated. Processing 
continues as normal. 

Now, in step 54, calculate Desired Energy by tools and reticles. Then, in step 56, fill in the 
missing values. This step is similar to the calculation of the tables in Phase I. However, Phase II is 
dynamic. New tools and new products will be released to production; thus, missing cells will be 
found and must be recalculated. 

In step 58, the Actual Energy is calculated. The correction value is calculated in step 60. 
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The APC system of the present invention has three modes: rapid mode, gradual mode, and 
generalized SPC mode. This is illustrated in Fig. 5. Gradual mode is designed to offset gradual 
drifts in the process while rapid mode is designed to do the same for sudden shifts in the process. 
The generalized SPC mode will choose between the two modes based on the deviation of the 
process output from the model predictions. That is, if the difference between the desired value and 
the actual value of the parameter is less than about 1-2 sigma in historical data calculation, the drift 
is said to be gradual. If the difference is more than about 3 sigma in historical data calaculation, the 
drift is said to be sudden. 

When running the first lot for a new mask, the APC system of the present invention will 
estimate energy by assigning the energy of the reticle which comes from the information 
technology (IT) system. Since this energy value is only an estimate, this lot must be excluded from 
the database. During running of the APC system of the present invention, four situations will 
result in a report: 

1. The process shift quantity is above a pre-determined value. 

2. The Desired Energy divided by the actual Energy is greater than a pre-determined value. 

3. An out of control situation occurs in the SPC monitor. 

4. An engineer desires to monitor the process. 

The APC process of the present invention has been implemented. With the process of the 
invention, the output CD is closer to the process target. More specifically, the rework rate was 
reduced by 80%. Throughput increased 0.43%, photoresist material savings were 0.31%, and 
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working hours were reduced by 15%. 

While the invention has been particularly shown and described with reference to the 
preferred embodiments thereof, it will be understood by those skilled in the art that various changes 
in form and details may be made without departing from the spirit and scope of the invention. 

What is claimed is: 
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